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Abstract: Using the gauge/string duality, we derive a set of Langevin equations describ- 
ing the dynamics of a relativistic heavy quark moving with constant average speed through 
the strongly-coupled AA = 4 SYM plasma at finite temperature. We show that the stochas- 
ticity arises at the string world-sheet horizon, and thus is causally disconnected from the 
black hole horizon in the space-time metric. This hints at the non-thermal nature of the 
fluctuations, as further supported by the fact that the noise term and the drag force in 
the Langevin equations do not obey the Einstein relation. We propose a physical picture 
for the dynamics of the heavy quark in which dissipation and fluctuations are interpreted 
as medium-induced radiation and the associated quantum-mechanical fluctuations. This 
picture provides the right parametric estimates for the drag force and the (longitudinal 
and transverse) momentum broadening coefficients. 
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1. Introduction 

Motivated by possible strong-coupling aspects in the dynamics of ultrarelativistic heavy 
ion collisions, there have been many recent applications of the AdS/CFT correspondence 
to the study of the response of a strongly coupled plasma — typically, that of the M = 4 
supersymmetric Yang-Mills (SYM) theory at finite temperature — to an external pertur- 
bation, so like a "hard probe" — say, a heavy quark, or an electromagnetic current (see the 
review papers |l|, [3| for details and more references). Most of these studies focused on 
the mean field dynamics responsible for dissipation (viscosity, energy loss, structure func- 
tions), as encoded in retarded response functions — typically, the 2-point Green's function 
of the N = 4 SYM operator perturbing the plasma. By comparison, the statistical prop- 
erties of the plasma (in or near thermal equilibrium) have been less investigated. Within 
the AdS/CFT framework, such investigations would require field quantization in a curved 
space-time — the AdS§ x S 5 Schwarzschild geometry dual to the strongly-coupled AA = 4 
SYM plasma — , which in general is a very difficult problem. Still, there has been some 
interesting progress in that sense, which refers to a comparatively simpler problem: that 
of the quantization of the small fluctuations of the Nambu-Goto string dual to a heavy 
quark immersed into the plasma. 

Several noticeable steps may be associated with this progress: In Ref. Q, a prescription 
was formulated for computing the Schwinger-Keldysh Green's functions at finite temper- 
ature within the AdS/CFT correspondence. With this prescription, the quantum thermal 
distributions are generated via analytic continuation across the horizon singularities in the 
Kruskal diagram for the AdS§ Schwarzschild space-time. Using this prescription, one has 
computed the diffusion coefficient of a non-relativistic heavy quark || , and the momentum 



- 1 - 



broadening for a relativistic heavy quark which propagates through the plasma at constant 
(average) speed || f?J. Very recently, in Refs. [||, a set of Langevin equations has been 
constructed which describes the Brownian motion of a non-relativistic heavy quark and of 
the attached Nambu-Goto string. Within these constructions, the origin of the 'noise' (the 
random force in the Langevin equations) in the supergravity calculations lies at the black 
hole horizon, as expected for thermal fluctuations. 

The Langevin equations in Refs. [| encompass previous results for the drag force 
|1C, 11] and the diffusion coefficient || of a non-relativistic heavy quark. But to our 
knowledge, no attempt has been made so far at deriving corresponding equations for a 
relativistic heavy quark, whose dual description is a trailing string flO| , Q |. In fact, the 
suitability of the Langevin description for the stochastic trailing string was even challenged 
by the observation that the respective expressions for the drag force and the momentum 
broadening do not to obey the Einstein relation J?]]. The latter is a hallmark of thermal 
equilibrium and must be satisfied by any Langevin equation describing thermalization. 
However, Langevin dynamics is more general than thermalization, and as a matter of facts 
it does apply to the stochastic trailing string, as we will demonstrate in this paper. 

Specifically, our objective in what follows is twofold: (i) to show how the Langevin 
description of the stochastic trailing string unambiguously emerges from the underlying 
AdS/CFT formalism, and (ii) to clarify the physical interpretation of the associated 
noise term, in particular, its non-thermal nature. 

Our main conclusion is that the stochastic dynamics of the relativistic quark is fun- 
damentally different from the Brownian motion of a non-relativistic quark subjected to 
a thermal noise. Within the supergravity calculation, this difference manifests itself via 
the emergence of an event horizon on the string world-sheet || (jj, which lies in between 
the Minkowski boundary and the black hole horizon, and which governs the stochastic 
dynamics of the fast moving quark. With our choice for the radial coordinate z in AdSs, 
the Minkowski boundary lies at z = 0, the black hole horizon at zjj = 1/T, and the world- 
sheet horizon at z s = Zh/^/j, where 7 = — v 2 is the Lorentz factor of the heavy 
quark. (We assume that the quark is pulled by an external force in such a way that its 
average velocity remains constant.) The presence of the world-sheet horizon means that 
the dynamics of the upper part of the string at z < z s (including the heavy quark at z ~ 0) 
is causally disconnected from that of its lower part at z s < z < zjj, and thus cannot be 
influenced by thermal fluctuations originating at the black hole horizon. 

This conclusion is supported by the previous calculations of the momentum broad- 
ening for the heavy quark H, [/J, which show that the relevant correlations are generated 
(via analytic continuation in the Kruskal plane) at the world-sheet horizon, and not at 
the black hole one. Formally, these correlations look as being thermal (they involve the 
Bose-Einstein distribution), but with an effective temperature T e & = T/^/7, which is the 
Hawking temperature associated to the world-sheet horizon. Thus, no surprisingly, our 
explicit construction of the Langevin equations will reveal that the corresponding noise 
terms arise from this world-sheet horizon. 

The Langevin equations for the relativistic heavy quark will be constructed in two 
different ways: (1) by integrating out the quantum fluctuations of the upper part of the 
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string, from the world-sheet horizon up to the boundary, and (2) by integrating out the 
string fluctuations only within an infinitesimal strip in z, from the world-sheet horizon at 
z = z s up to the 'stretched' horizon at z = z s (l — e) with e <C 1; this generates a 'bulk' 
noise term at the stretched horizon, whose effects then propagate upwards the string, via 
the corresponding classical solutions. Both procedures provide exactly the same set of 



Langevin equations, which encompass the previous results for the drag force |1(], 11] and 
for the (longitudinal and transverse) momentum broadening |(| |?]]. In these manipulations, 
the lower part of the string at z > z s and, in particular, the black hole horizon, do not play 
any role, as expected from the previous argument on causality. 

If the relevant fluctuations are not of thermal nature, then why do they look as being 
thermal ? What is their actual physical origin ? And what is the role played by the 
thermal bath ? To try and answer such questions, we will rely on a physical picture for 
the interactions between an energetic parton and the strongly-coupled plasma which was 
proposed in Refs. [12, 13, [14], |2|, and that we shall here more specifically develop for the 
problem at hand. In this picture, both the energy loss ('drag force') and the momentum 
broadening ('noise term') are due to medium-induced radiation. This is reminiscent of 
the mechanism of energy loss of a heavy, or light, quark at weak coupling [15, 16, 17, 18], 
with the main difference being in the cause of the medium-induced radiation. At weak 
coupling, multiple scattering off the plasma constituents frees gluonic fluctuations in the 
quark wavefunction, while at strong coupling the plasma exerts a force, proportional to T 2 , 
acting to free quanta from the heavy quark as radiation. In the gravity description, this 
appears as a force pulling energy in the trailing string towards the horizon. At either weak 
or strong coupling, quanta are freed when their virtuality is smaller than a critical value, 
the saturation momentum Q s ; at strong coupling and for a fast moving quark, this scales 
like Q s ~ \f^T ■ Within this picture, the world-sheet horizon at z s ~ 1/Q S corresponds 
to the causal separation between the highly virtual quanta (Q ^> Q s ), which cannot decay 
into the plasma and thus are a part of the heavy quark wavefunction, and the low virtuality 
ones, with Q < Q s , which have already been freed, thus causing energy loss. The recoil of 
the heavy quark due to the random emission of quanta with Q < Q s is then responsible 
for its momentum broadening. 

From his perspective, the noise terms in the Langevin equations for the fast moving 
quark reflect quantum fluctuations in the emission process. Of course, the presence of 
the surrounding plasma is essential for this emission to be possible in the first place (a 
heavy quark moving at constant speed through the vacuum could not radiate), but the 
plasma acts merely as a background field, which acts towards reducing the virtuality of the 
emitted quanta and thus allows them to decay. The genuine thermal fluctuations on the 
plasma are unimportant when 7 > 1, although when 7 ~ 1 they are certainly the main 
source of stochasticity, as shown in [g, [|. Besides, we see no role for Hawking radiation of 
supergravity quanta at any value of 7. 

This picture is further corroborated by the study of a different physical problem, where 
the thermal effects are obviously absent, yet the mathematical treatment within AdS/CFT 
is very similar to that for the problem at hand: this is the problem of a heavy quark 
propagating with constant acceleration a through the vacuum of the strongly-coupled 
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N = 4 SYM theory fl4|, 16, p(T]. The accelerated particle can radiate, and this radiation 



manifests itself through the emergence of a world-sheet horizon, leading to dissipation and 
momentum broadening. The fluctuations generated at this horizon are once again thermally 
distributed, with an effective temperature T e g- = a/2n. In that context, it is natural to 



interpret the induced horizon as the AdS dual of the Unruh effect [21] : the accelerated 



observer perceives the Minkowski vacuum as a thermal state with temperature a/2n. For 
an inertial observer, this is interpreted as follows p^] : the accelerated particle can radiate 
and the correlations induced by the backreaction to this radiation are such that the excited 
states of the emitted particle are populated according to a thermal distribution. Most likely, 
a similar interpretation holds also for the thermal-like correlations generated at the world- 
sheet horizon in the problem at hand — that of a relativistic quark propagating at constant 
speed through a thermal bath. It would be interesting to identify similar features in other 
problems which exhibit accelerated motion, or medium-induced radiation, or both, so like 
the rotating string problem considered in Ref . 22 ] . 

The paper is organized as follows: In Sect. 2 we construct the Langevin equations de- 
scribing the stochastic dynamics of the string endpoint on the boundary of AdS§, i.e., of the 
relativistic heavy quark. Our key observation is that, in the Kruskal-Keldysh quantization 
of the small fluctuations of the trailing string, the stochasticity is generated exclusively 
at the world-sheet horizon. This conclusion is further substantiated by the analysis in 
Sect. 3 where we follow the progression of the fluctuations along the string, from the 
world-sheet horizon up to the string endpoint on the boundary. We thus demonstrate that 
the noise correlations are faithfully transmitted from the stretched horizon to the heavy 
quark, via the fluctuations of the string. Finally, Sect. 4 contains our physical discussion. 
First, in Sect. 4.1, we argue that the Langevin equations do not describe thermalization, 
although they do generate thermal-like momentum distributions, but at a fictitious tem- 
perature which is not the same as the temperature of the plasma, and is moreover different 
for longitudinal and transverse fluctuations. Then, in Sect. 4.2, we develop our phys- 
ical picture for medium-induced radiation and parton branching, which emphasizes the 
quantum-mechanical nature of the stochasticity. 



2. Boundary picture of the stochastic motion 

In this section we will construct a set of Langevin equations for the stochastic dynamics 
of a relativistic heavy quark which propagates with uniform average velocity through a 
strongly-coupled J\f = 4 SYM plasma at temperature T. To that aim, we will follow the 
general strategy in Ref. Q, that we will extend to a fast moving quark and the associated 
trailing string. In this procedure, we will also rely on previous results in the literature 
H 0] concerning the classical solutions for the fluctuations of the trailing string and their 
quantization via analytic continuation in the Kruskal plane. 

2.1 The trailing string and its small fluctuations 

The AdS dual of the heavy quark is a string hanging down in the radial direction of 
AdSs, with an endpoint (representing the heavy quark) attached to a D7-brane whose 
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radial coordinate fixes the bare mass of the quark. The string dynamics is encoded in the 
Nambu-Goto action, 



S = J d 2 cr v / -det h a/3 , h al 3 = 9^9^ d^x" , (2.1) 

where a a , a = 1, 2, are coordinates on the string world-sheet, h a p is the induced world- 
sheet metric, and g^v is the metric of the ^cfSs-Schwarzschild space-time, chosen as 

R 2 ( ., , , 2 dz 2 



d ' =^{-^ dt +ix + W))' (2 - 2) 

where f(z) = 1 — z 4 and T = 1/ttzh is the Hawking temperature. (As compared to the 
Introduction, we have switched to a dimensionless radial coordinate.) 

The quark is moving along the longitudinal axis x 3 with constant (average) velocity 
v in the plasma rest frame. For this to be possible, the quark must be subjected to some 
external force, which compensates for the energy loss towards the plasma. The profile of 
the string corresponding to this steady (average) motion is known as the 'trailing string'. 
This is obtained by solving the equations of motion derived from Eq. fl2.l| ) with appropriate 



boundary conditions, and reads [Uj, 11] 



Xq = v t + — ( arctan z — arctanhz) . (2-3) 

In what follows we shall be interested in small fluctuations around this steady solution, 
which can be either longitudinal or transverse: x s = Xq + Sxe(t, z) and x± = 5x±(t, z). To 
quadratic order in the fluctuations, the Nambu-Goto action is then expanded as (in the 
static gauge a a = (t, z)) 

S = - ^J** f dtdz\ + f dtdzP a d a 5x e 



- ~ J dtdz Tf d a 5x l d p 5x l + Tf d a 5x\d p 5x\ , (2.4) 
where z s = \/\ — v 2 = l/y^y and 1 @ 

P^ = H^1(^)) : (2.5) 



T 



.a 3 



ArpClf} 



TT 



Vat 2 



2z 2 



Z H l ~( ZZ ->) A V 2 

~z? (l-z*) 2 T^z* 



(2.6) 



1-z 4 



The quantities T"^ have the meaning of local stress tensors on the string. At high energy, 
the components of the longitudinal stress tensor are parametrically larger, by a factor 



1 In Eq. (2.6), we have corrected an overall sign error in Eq. (22) of Ref. 
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7 2 3> 1, than the corresponding components of the transverse stress tensor. This 
difference reflects the strong energy-dependence of the gravitational interactions. 

Using d a P a = 0, one sees that the term linear in the fluctuations in Eq. ( |2.4| ) does not 
affect the equations of motion, which therefore read 

d a (T#dffll,) = 0, V = 8x e , Sx ± , (2.7) 

in compact notations which treat on the same footing the longitudinal and transverse 
fluctuations. Upon expanding in Fourier modes, 

^u,z)e-^\ (2.8) 

-oo 27T 

this yields 

{a{z)d 2 z - 2b(uj, z)d z + c(w, z)} ip(uj, z) = 0, (2.9) 

where 

a(z) = z(l-z 4 f{z 4 -z% 
b(u, z) = (1 - z 4 ) [l - z 8 - v 2 (l - z 2 + iuJz H z 3 )] , 
c(uj, z) = ujzhz \u)Zh(1 — z 4 ) + v 2 z 4 {ujzh + 4iz)] ■ (2-10) 

The zeroes of a(z) determine the regular singular points of this equation. In particular, 
the special role played by the point z s as a world-sheet horizon becomes manifest at this 
level: for z = z s the value of d z ip(uj,z s ) is determined from the equation of motion. This 
means that fluctuations of the string at z < z s are causally disconnected from those below 
the location of the world-sheet horizon. 

2.2 Keldysh Green function in AdS/CFT 

In what follows we construct solutions to Eqs. (|2.9| )- (|2.10| ) for the string fluctuations which 
are well defined everywhere in the Kruskal diagram for the AdS§ Schwarzschild space-time 
(see Fig. [j]). These solutions are uniquely determined by their boundary conditions at the 
two Minkowski boundaries — in the right (R) and, respectively, left (L) quadrants of the 
Kruskal diagram — , together with the appropriate conditions of analyticity in the Kruskal 
variables U and V (as explained in Q). The latter amount to quantization prescriptions 
which impose infalling conditions on the positive-frequency modes and outgoing conditions 
on the negative-frequency ones. These prescriptions ultimately generate the quantum 
Green's functions at finite-temperature and in real time, which are time-ordered along the 
Keldysh contour Specifically, the time variables on the R and, respectively, L boundary 
in Fig. H correspond to the chronological and, respectively, antichronological branches of 
the Keldysh time contour. 

As usual in the framework of AdS/CFT, we are interested in the classical action ex- 
pressed as a functional of the fields on the boundary. We denote by ipu(tfi, z) and i^l^Ri z) 
the classical solutions in the R and L quadrant, respectively. Making use of the equations 
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Figure 1: Kruskal diagram for AdS$ Schwarzschild metric; the position of the induced horizon on 
the string world-sheet is shown with dashed lines in both R and L quadrants. 



of motion and integrating by parts, the classical action reduces to its value on the boundary 
of the Kruskal diagram, i.e., the R and L Minkowski boundaries: 



Si 



bndry 



dt R 



ZR=Zm 



dt L 



(2.11) 



where it is understood that the terms involving P z exist only in the longitudinal sector. 
z m <C 1 is the radial location of the D7-brane on which the string ends. 



In Eq. ( 2.11 ), the world-sheet index (3 can take a priori both values t and z, but the 
contribution corresponding to (3 = t is in fact zero, since the respective integrand is an 
odd function of t. This is worth noticing since in Ref. ^ it was found that the dominant 
contribution to the imaginary part of the retarded propagator at high energy (7 > 1) 
comes from the piece proportional to T zt . We will later see how that contribution arises 
in the present calculation, where only the piece proportional to T zz survives in Eq. ( 2.1 1| ) . 

Switching to the frequency representation, we introduce a basis of retarded and ad- 
vanced solutions, ip ret (uj,z) and i() a( i v (uj, z), which are normalized such that ip re t(uj,0) = 
^advi^, 0) = 1. They obey ip re t(uj, z) = ip* et (—u>, z), and similarly for tp a dv These solutions 
are truly boundary-to-bulk propagators in Fourier space. They have been constructed in 
Ref. Q (see also @) from which we quote the relevant results. 

Note first that, unlike what happens for a static (or non-relativistic) quark |J, the re- 
tarded and advanced solutions are not simply related to each other by complex conjugation: 
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one rather has 



iW", *) = bwr /2 b(^/^)]-^ /22s w* e *K , (2.i2) 

with g(z) = j±| e -2arctan2_ N ear the boundary [z -C 1), these solutions behave as follows 



VVet(w,z) 



2 2 \ 
1 + ^"Z 2 + 0{*)\ + C ret (oo) {z 3 + 0(Z 5 )) (2.13) 



^«fo(">*) = (l + + 0(^ 4 )) + C^MO* 3 + 0(z 5 )), (2.14) 

where the expansion involving even (odd) powers of z is that of the non-normalizable 
(normalizable) mode. The coefficients of the normalizable mode are related by 

C adv (u) = C? et (u)-iX(uj), (2.15) 

with 

X(u) = ^ vW , ImC re( ( W ) = i ^ . (2.16) 



The real part of coefficient C re t{oj) has been numerically evaluated in Ref. [jTj. Here, we 
only need to know that, at small frequency uj <C z s T, its real part is comparatively smaller: 
ReC re f(w) ~ O(L0 2 zfj/Zg). Note that, at high energy (7 1), X dominates over lm.C re t 
in Eq. <&T$ ). 

Consider also the approach towards the world-sheet horizon [z = z s ) from the above 
(z < z s ): there, tp re t remains regular, whereas VWu has a branching point: 

iPadv{u,z) (X (z s - z)^ [1 + 0(z s - z)} . (2.17) 

This shows that ^ a d v (uJ,t, z) = e~~ lujt ip a d v (u>, z) is an outgoing wave: with increasing time, 
the phase remains constant while departing from the horizon. One can similarly argue that 
^S> re t(uj,t, z) = e~ lujt ipret(u, z) is an infalling solution |@]. 

We now expand the general solution in the right and left quadrants of the Kruskal 
diagram in this retarded/advanced basis : 

1pR(u,z) = A(Lv)lpret(uJ,z) + B(ui)t/j adv (u, z) , 

ip L {u],z) = C{ui)ip ret {uJ,z) + D(u)ip adv (uj,z). (2.18) 

We need four conditions to determine the four unknown coefficients A, B, C, and D. Two 
of them are provided by the boundary values at the R and L Minkowski boundaries, that 
we denote as and ip^tv), respectively. The other two are determined by analyticity 

conditions in the Kruskal plane, which allows one to connect the solution in the L quadrant 
to that in the it! quadrant. With reference to Fig. |], one sees that this requires crossing 
two types of horizons: the world-sheet horizons in both R and L quadrants, and the 
black hole horizons at XJ = and V = 0. The detailed matching at these horizons, 
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following the prescription of Ref . || , has been performed in Appendix B of Ref . Q , with 
the following results: the multiplicative factors associated with crossing the black hole 
horizons precisely compensate each other, unlike those associated with crossing the world- 
sheet horizons, which rather enhance each other. Hence, the net result comes from the 
world-sheet horizons alone 2 , and reads 

/ i n \ / a \ 

(2.19) 

The two independent coefficients can now be determined from the boundary values %I) r {uj) 
and tp^(io). This eventually yields 

A(oj) = (1 + n)(cu)V^M " n(oj)^° L (u), (2.20) 





B(u) = n(u) - V^M] . (2.21) 

Here n{uj) = \/(e u) l ZsT — 1) is the Bose-Einstein thermal distribution with the effective 
temperature T e g = z s T = T / ^7. As it should be clear from the previous manipulations, 
this effective thermal distribution has been generated via the matching conditions at the 
R and L world-sheet horizons, cf. Eq. ( |2.19[ ). 

The equations simplify if one introduces 'average' (or 'classical') and 'fluctuating' vari- 
ables, according to ip r = (ip R + V'l)/2 and ip a = ip R — ipL, and similarly for the boundary 
values. One then finds 

1p r (uJ, z) = 1p°(u)lp ret (uJ, z) + 1 + ^L^_ -^(^(^(o;, z ) - l{) adv (uj, z)), 

^a(w,s) = (2.22) 
and the boundary action takes a particularly simple form: 

Sbndry = ^J T zz \z)[^ r {-u , z)d z ^ a {uj , z) + (r o)]^. (2.23) 

(We have here omitted the term linear in the fluctuations, since this does not matter for 
the calculation of the 2-point Green's functions. This term will be reinserted in the next 
subsection.) By combining the above equations, we finally deduce 

^bndxy = ~i J ^ $(-a;)G&(a;)$( W ) ~\\% V> a (-^)G sym MV», (2.24) 
with the retarded and symmetric Green's functions defined as 

G^ R {u) = z\G\ R (u) = -\Ti z (z)d z U ret (u,z)ij> adv (-u,z)\ , (2.25) 



2 This point is even more explicit in the analysis in Ref. M, where a different set of coordinates was used, 
in which the world-sheet metric is diagonal. With those coordinates, the only horizons to be crossed when 
going from the R to the L boundary in the respective Kruskal diagram are the world-sheet horizons. 
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and, respectively, 

G sym (u;) = -(l + 2n(cj))ImG° R (uj). (2.26) 
Note that Eq. ( p^26| ) is formally the same as the fluctuation-dissipation theorem (or 'KMS 



condition') characteristic of thermal equilibrium, but with an effective temperature T c g = 
T. By also using Eqs. ( |2.13| )-( |2.14 ) together with the expression of TJ 2 given in (p.6|), 



Zs 

one finally deduces 



G±,hM = = G r (oj) - jMqu" , (2.27) 

where 

G R (u) = -Y[C ret {u) + C* adv {u)), 

3vA . , vAT -v/A^m 

y = 5 , M Q = = £ . (2.28) 

Am^z\ w 2z m 2ttR 2 v ' 

Mq is the (bare) rest mass of the heavy quark and is independent of temperature, as 
manifest in his last rewriting. (r m denotes the position of the D7-brane in the usual radial 
coordinate r, which is related to z as z/irT = R 2 /r.) At finite temperature, this mass 
receives thermal corrections, as encoded in the contribution of 0{uj 2 ) to KeC re t(uj); such 
corrections are however negligible at high energy, since their contribution to G r (lo) is not 
enhanced by a factor of 7 (unlike Mq). 

Note that the previous formulae fully specify the imaginary part of the retarded prop- 
agator, and hence also G sym (u;). Namely, by using (cf. Eqs. ( 2.15| )-( p.l6| )) 



ImCrM + lmC:^) = ^(1 + ,V) = ^7 2 , (2-29) 
one immediately finds 

lmG R (u}) = -ujrj, with rj = — y- T 2 . (2.30) 

Remarkably, this exact result involves just a term linear in oj. On the other hand, we 
expect KeG R (uj) to receive contributions to all orders in uj starting at 0(lo 2 ). 

The above expression for G R , Eq. ( |2.27| ), coincides with that originally derived in 
Ref. @, although the definition used there for the retarded propagator was different, namely 

G° R (u) = -% et T^d^ ret \ z=Zm . (2.31) 

(Recall that ^> re t{uo, t, z) = e~ tujt ip re t(uj, z).) With this definition, the dominant contri- 
bution to the imaginary part at high energy — the term proportional to X(uu) — arises 
from the time derivative of the retarded solution. Although it does not naturally emerge 



when constructing the boundary action in the Kruskal plane, this formula ( 2.31 ) has an- 
other virtue, which will be useful later on: with this definition, the imaginary part of the 
retarded propagator, 

ImG^H = j,T^[% et dp^ ret - *r*dp*U) . ( 2 - 32 ) 
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can be evaluated at any z, since the r.h.s. of Eq. ( 2.32 ) is independent of z. Indeed, as 
noticed in Ref. |7]], the world-sheet current 

ja = Yi T ^(^*soid^soi - *soldf3** sol ) , (2.33) 



(^> so i(t, z) is an arbitrary solution to the classical EOM (|2.7| )) is conserved by the equations 
of motion, d a J a = 0. When ^! so i(t,z) = ^ re t(uj,t, z), this conservation law reduces to 
d z J z = 0. As we shall shortly see, lmG R (uj) is a measure of the energy loss of the heavy 
quark towards the plasma. Thus the fact this quantity is independent of z is a statement 
about the conservation of the energy flux down the string in the present, steady, situation. 

2.3 A Langevin equation for the heavy quark 



Following the general strategy of AdS/CFT, the boundary action ( p. 24 ) can be used to 



generate the correlation functions of the M = 4 SYM operator which couples to the bound- 
ary value of the field — in this case, the Schwinger-Keldysh 2-point functions of the force 
operator acting on the heavy quark 0, [6|, [5], 0]. Alternatively, in what follows, this action 
will be used to derive stochastic equations for the string endpoint, in the spirit of the 
Feynman- Vernon 'influence functional' [^4|] (see also Ref. ||). 

To that aim we start with the following path integral which encodes the (quantum and 
thermal) dynamics of the string fluctuations in the Gaussian approximation of interest: 

Z = j [DiP° R D^ R ] [DiP° L D?p L ] e iS *- iS L . (2.34) 

This involves two types of functional integrations: (i) those with measure [D^rDi^l], 
which run over all the string configurations ipR,L(t,z) (in the corresponding quadrants of 
the Kruskal diagram) with given boundary values V^lW' an d (ii) those with measure 
[Dip R Dip^], which run over all the possible paths iPrl^) for these endpoint values. 

Performing the Gaussian path integral over the bulk configurations amounts to eval- 
uating the action in the exponent of Eq. ( 2.34[) with the classical solutions computed in 



the previous section. This leaves us with the boundary action in Eq. ( 2.24j ), which deter- 



mines the dynamics of the string endpoints — i.e., of the heavy quark — , and which is 
itself Gaussian. To perform the corresponding path integral it is convenient to 'break' the 
quadratic term for the fluctuating fields ip®i by introducing an auxiliary stochastic field 
£(i). Then the partition function becomes 



exp Li [ dMi>°(t)[G°e(M')Vv (O + <*(* " 0(** " £(0)] \ » ( 2 -35) 



where we recall that the term involving P z appears only in the longitudinal sector. The 
integral over ip a acts as a constraint which enforces a Langevin equation for the 'average' 
field ip r . This equation reads 

J dt'G° R (t, 0$(0 +P Z ~ ^) = 0, <£(t)£(0> = G sym(t, , (2.36) 
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and is generally non-local in time. At this point it is convenient to focus on the large time 
behaviour, as controlled by the small-frequency expansion of the Green's functions Gr and 
G sym , and also distinguish between longitudinal and transverse fluctuations. As discussed 
in Sect. 2.2, for oj <C z s T, the retarded propagator reduces to its imaginary part, Eq. ( pE3C| ) 
(in addition to the bare mass term). In the same limit, one can use 1 + 2n(u) ~ 2z s T/u) 
to simplify the expression of G sym (tL>), which then becomes independent of u : 

G ± , sym (uj) ~ ttVXj^T 3 = K ± , 

G i>sym (uj) ~ ttVXj^T 3 = Kt . (2.37) 

This in turn implies that, when probed over large time separations t — t' 3> l/z s T, the 
retarded propagator can be replaced by a local time derivative ('friction force'), while the 
noise-noise correlator looks local in time ('white noise'). The we can write 



d 2 Sx± dSx 



iMq ^d? = ~ ir] ~^it + <^(*)^-l(*')> = - *o > ( 2 - 38 ) 

for the transverse modes and, respectively (note that P z = jrjv), 

7 3 m q ^ = - 7 3 v ^ - -ym + m, &(t)6(0> = s(t - 1') , (2.39) 

for the longitudinal one. The physical interpretation of these equations becomes more 
transparent if they are first rewritten in terms of the respective momenta p± = jMqv± 
and pi = jMqvi, with v± = dSx±/dt and v i = v + dSxn/dt. 

At this point, we come across a rather subtle point: in all the equations written so far, 
the Lorentz factor 7 is evaluated with the average velocity v of the heavy quark — the one 
which enters the trailing string solution (2.3). However, the event-by-event fluctuations 
of the velocity turn out to be significantly large (especially in the longitudinal sector; see 
below), and then it becomes appropriate to define the event-by-event (or 'fluctuating') 
momenta p± and p£ by using the respective, event-by-event, Lorentz factor, as evaluated 
with the instantaneous velocity. For more clarity, let us temporarily denote by and 70 
the average velocity and the associated Lorentz factor, 70 = l/y'T— Vq, and reserve the 
notations v and 7 for the respective fluctuating quantities: 

When taking the time derivatives of p± and p#, as associated with variations in v± and, 
respectively, vg, one must also take into account the corresponding change in the 7-f actor. 
Consider the longitudinal sector first: 



1 dpi 
Mr, dt 



^)^(^V)^^, (2.4!) 

where the last, approximate, equality follows since the fluctuations are assumed to be small, 
hence ve — vq and 7 ~ 70. The final result above is recognized as the expression in the 
l.h.s. of Eq. ( 2.39| ). To the same accuracy, we can write (with 5vg = d5xn/dt) 



7Vl ~ ( 7o + — Sve ) (vo + Sv e ) ~ 7 u + (70 + 7o v o)<^ = 7o^o + lo Sv e > ( 2 - 42 ) 
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in which we recognize the terms multiplying 77 in the r.h.s. of Eq. flp|). 
Consider similarly the transverse sector. The analog of Eq. fl2,41| ) reads 

1 dpx / \dv ± f , . 2 2\ d 2 fe_L d 2 <5x_L , . 

^ + v ^7^)^T = ^( 1 + v ^)^u^^^-n 2 -, (2.43) 



M Q dt V dvj_J dt ' " Xl > dt 2 - ,u dV ' 

where we made the additional assumption that v\ <C 1 — Vq. (This can be always ensured 
by taking the quark mass Mq to be sufficiently large.) Similarly, in the r.h.s. of Eq. ( 2.3§| ), 
one can replace 7oVj_ ~ P±/Mq. 

To summarize, to the accuracy of interest, we have derived the following Langevin 
equations for the dynamics of the heavy quark 

d// ' -VDPl + £(t), (fi(t)ei(f )) = K±^S(t - 1') , (2.44) 



dt 

-VDPe + te(t), Mt)Ut')) = nt 5(t - t') , (2.45) 



dpi 



dt 



where the upper index i = 1, 2 in Eq. ( 2.44| ) distinguishes between the two possible trans- 



verse directions, kj_ and Kg are given in Eq. (EE37D , and 

_ t] vr\/A 2 
^ = W Q = 2M- Q T - (2 - 46) 

The general structure of these equations — with a friction term (or 'drag force') describing 
dissipation and a noise term describing momentum broadening — is as expected, and so 
are the above expressions for 77/5, k± and K£, which agree with previous calculations in 
the literature [||, ||, [jj. It is however important to keep in mind that Eqs. (2.44)-( 2.46| ) 



have been derived here only for the situation where the fluctuations in the velocity of the 
heavy quark remain small as compared to its average velocity vq. To ensure that this 
is indeed the case, Eq. ( 2.45] ) for the longitudinal motion must be supplemented with an 



external force which is tuned to reproduce the average motion. (Without such a term, 
Eq. ( 2.45| ) would describe the rapid deceleration of the heavy quark due to its interactions 



in the plasma. Such a deceleration may entail additional phenomena, like bremsstrahlung, 



which are not encoded in the above equations; see the discussion in Refs. 14, [L9|, |23j1 .) 



Namely, we shall add to the r.h.s. of Eq. ( [2,45 ) a term F ex t = i]^qVo which for large times 



equilibrates the average friction force and thus enforces a constant average velocity vq. 
Further consequences of these equations will be discussed in Sect. 4. 

3. Bulk picture of the stochastic motion 

In the previous section we have obtained a set of Langevin equations for the heavy quark 
by integrating out the fluctuations of the upper part of the string, from the world-sheet 
horizon up to the boundary. The noise terms in these equations have been generated via 
boundary conditions at the world-sheet horizon, cf. Eq. ( p. 19 ). This suggests that, within 



the context of the supergravity calculation, quantum fluctuations are somehow encoded 
in the world-sheet horizon. To make this more explicit, we shall follow Refs. ]9|, [8| and 
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construct a set of equations describing the stochastic dynamics of the upper part of the 
string, in which the noise term is acting on the lower endpoint, infinitesimally close to the 
world-sheet horizon. 

More precisely, we introduce a 'stretched' horizon at Zh = z s — e and integrate out the 
fluctuations of the part of string lying between z s and z^. The procedure is quite similar to 
the one described in the previous section except that one has to fix the boundary values for 
the fluctuations also on the stretched horizon, rather than just on the Minkowski boundary. 
Denoting the respective values by ip h , where as before ip stands generically for either 5x£ 
or 5x± , this procedure yields an effective action S^ s for tp h with the same formal structure 



as exhibited in Eq. ( 2.24 ), that is, 



iS " s = ~ i j |^(-^MVV>) -\J ^(-^G^ico)^)- (3-1) 

(We temporarily omit the term linear in the fluctuations; this will be restored in the final 
equations.) The horizon Green's functions G R and Gg ym will be shortly constructed. The 
calculations being quite involved, it is convenient to start with a brief summary of our main 
results: 

The r-fields tp r (cj,z) describing the string fluctuations within the bulk [z m < z < 
Zh) obey the equations of motion (2/7) with Neumann boundary conditions at z = z m 



- meaning that the string endpoint on the boundary is freely moving (except for the 
imposed longitudinal motion with velocity vq) — and with Dirichlet boundary conditions 
at z = z^ if) r (u,Zh) = ipr(u>). This boundary field ipri^) is however a stochastic variable, 
whose dynamics is described by the effective action ( |3.1| ). Via the classical solutions, this 
stochasticity is transmitted to the upper endpoint of the string, i.e., to the heavy quark. 
As a result, the latter obeys the same Langevin equations as previously derived in Sect. 2. 

We start with the partition function encoding the quantum dynamics of the upper part 
of the string (z m < z < Zh) in the Gaussian approximation: 

Z = j D^ R Di> R D4) h R Di)° L D^ L D^ h L e i5 «- i5i+i5 eff. (3.2) 

The different measures D^°, Dip and Dip h correspond, respectively, to the path integral 
over the string endpoint on the Minkowski boundary, over the bulk part of the string, 
and over the string endpoint on the stretched horizon (separately for the left and right 



quadrants of the Kruskal plane). Also, Sr and Sl are defined as in Eq. (2.4), but with the 
integral over z restricted to z m < z < z^. In what follows we shall construct the various 
pieces of the action which enter the exponent in Eq. (|3.2| ) . 

(I) The effective action at the stretched horizon, S% s . As anticipated, this is obtained by 
integrating out the string fluctuations within the infinitesimal strip Zh < z < z s . To that 
aim, we need the classical solutions ipn(uj,z) and x/jl(uj,z) in the Kruskal plane which take 
the boundary values ip R {u>) and V'l( u; ) a ^ 2 = z h an d are related by the condition ( 2.1S| ). 
Clearly, the respective solutions read (in the (r, a) basis, for convenience) 

T/V(W, z) = WV)^retK Z) + 1 + VaC^X^Uw, *) " V&faO", *)). 

1> a (u>,z) = ^(u)^ v (u,z), (3.3) 



-14- 



where ip^ et and ip^dv are resca l e d versions of the retarded and advanced solutions introduced 
in Sect. 2.2 which are normalized to 1 at z = Zh] e.g., ^ et (u},z) = ip re t(uj, z) / ip re t(ui , Zh). 
For z close to z a (and hence to Zh as well), these functions can be expanded as 



1 + 0(z s 



z), 

2z s 



[l + 0{z s -z)\ . 



(3.4) 



Substituting these classical solutions into the action produces the boundary action shown 
in Eq. ([O]) , with G R defined by the horizon version of Eq. ( |2.25[ ) . Given the near-horizon 
behaviour of the solutions (3.4) and of the local tension T zz (z) (which vanishes at z = z s , 
cf. Eq. fl2.6| )), it is clear that only d z tp^ dv contributes to G R in the limit e — » 0. This yields 
the following, purely imaginary, result: 



G h ± , R (u) = z A s Gl R (u) = ~TT(z h )d z ^(-u,z) 



z=z h 



(3.5) 



This coincides, as it should, with the imaginary part of the respective boundary propaga- 
tor 3 , Eq. ( 2.30P (cf. the discussion at the end of Sect. 2.2). Then the symmetric Green's 
function G^ ym , which is related to ImG^ via the KMS relation ( |2.26| ), is exactly the same 
as the corresponding function on the boundary. 

If the string point on the stretched horizon was a free endpoint endowed with the action 
( |3.1| ), it would obey Langevin equations similar to those derived in Sect. 2.3. However, this 
is an internal point on the string, and as such it is also subjected to a tension force from 
the upper side of the string at z < Zh- This force is encoded in the bulk action Sr — Sl, 
to which we now turn. 

(II) The bulk piece of the action Sr — Sl- This is defined as 



Sr — Sl 



dz / dtrt 



■-h 



rpRd a (T^d^-rp L d a (T a ^dpi; L 
dtT zz (z)U R d z <<p R -i; L d z ip L 



dz / dt 



z=z h 



(3.6) 



or, after going to Fourier space and to the (r, a)-basis, 



Sr — Sl 



doj 
2^ 



dzip a (-uj,z)d a T af3 (z)dpip r (u,z) 



2 / -^ TZZ ( z )(M-u,z)d z M",z) + (r «-> a) 



Z = Zh 



(3.7) 



We were so explicit here about the integration by parts, because this operation turns out 
to be quite subtle. First, notice that the contributions proportional to T zt have cancelled 



incidentally, this calculation of ImG^, which is exact, together with the conservation law d z J z = 0, 
cf. E q. (p.33| ), can be used to check, or even derive, the expressions for ImC re t(ui) and Im C a dv(uj) given in 



Eqs. (|2.1E|)-(^1^). 
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in the boundary terms, for the same reason as discussed below Eq. ( |2.11 ), i.e., because they 
are odd functions of t (or lo). To ensure this property, it has been important to perform 
the previous operations in the order indicated above, that is, to first integrate by parts, as 
in Eq. ( |3.6| ), and only then change to the (r, a)-basis, as in Eq. (p.7|). (Reversing this order 
would have affected the symmetry properties of the integrand, and then the terms oc T zt 
would not cancel anymore.) 

Second, there are some subtleties about the boundary value of Sr — Sl at the stretched 
horizon, that we rewrite here for more clarity: 



(Sr - Sx)b; 



bndry 



do; 

2^ 



(z h )(^tp a (-u;,z)d z ip r (u;,z) + (r 



(3.8) 



= 2/1 



If we were to evaluate this action with the classical solutions (3.3), the result would precisely 
cancel the effective action ([O]) in the exponent of Eq. d3.2[ ). Indeed, up to a sign, Eq. ( |3~^ ) 
has exactly the structure that has been used to build the effective action by inserting the 
classical solutions (compare to Eq. fl2.23Q ). However, in the present context, Eq. (|3.8|) must 
be rather seen as the boundary value of the bulk action when approaching the stretched 
horizon from the above (i.e., from z < z^), and as such it provides boundary conditions for 
the dynamics of the upper side of the string (see below). This being said, it is nevertheless 
possible, and also convenient, to use the appropriate piece of Eq. (|3.8|) in order to cancel 
the dissipative piece oc G R in the effective action ( |3.1[ ). This is simply the statement that 
the totality of the energy which crosses the stretched horizon coming from the above flows 
further down along the string. 

Specifically, the relevant piece of Eq. Q3.8|) is that proportional to d z ip a , which after 



using Eq. (3.3) can be evaluated as 

{z h )d z if)% dv (-uj,z 
where we have recognized the expression 



\_rj<ZZ I 

2 



=2;, 



5|) for G R . One thus obtains: 



(Sr - Sl)^ 



bndry 



duj 
2^ 



■T"(z)d z il> T {u,z) 



= 2/. 



As anticipated, the last term in Eq. ( |3.10|) compensates the piece involving G R in Eq. (| 
and then the total action reads 



(3.9) 

(3.10) 

), 



iS R - iS L + iS; 



eff 



du 



+ % 



2tt 



-u 



T zz (z) d z M^, z) + ^(u)d z ip adv (-cj, z) 



pz 



dz 



^Va(-^)d a \T^(z)d^ r (u,z) 



duJ tht n 
2^^ ( " W) 



rpZZ i 

2 



z)dMu, Z )-e(u) 



(3.11) 



= 2h 



where the differences between longitudinal and transverse fluctuations are kept implicit 
(in particular, it is understood that the term proportional to P z appears only in the 
longitudinal sector) . Two additional manipulations have been necessary to write the action 
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in its above form: (i) In the first line of Eq. ( |3.1lD we have used ^ a (w, z) = ^{oJ^advi^i z ), 
cf. Eq. ( 2,22 ). (ii) The piece involving G^ ym (w) in Eq. (3A) have been reexpressed as a 



Gaussian path integral over the noise variables £ , which therefore obey 



(£ h (u)£ h (u/)) = 2vr£(w + u/) (1 + 2^77/ , with n(u) 



1 



3 uj/z s T _ ^ 



(3.12) 



Once again, Eq. ( 3.12 ) involves the effective temperature T e fj = z s T. 

By integrating over the fluctuating fields ip®(—uj), ip a (—u,z), and a;), we are 
finally left with the following set of equations of motion and boundary conditions: 

(1) A modified Neumann boundary condition for the string endpoint string at the 
boundary (we temporarily reintroduce the polarization label p with p = t or _L) : 



-T; z {z) d z r r ^,z)+^(u;)d z ^ adv (-u;,z) 



T]JV d p i . 



(3.13) 



(2) The standard equations of motion for the fluctuations i/j r (uj, z) of the string in the 
bulk at z m < z < z h (cf. Eq. pTzj)). 

(3) A stochastic equation for the string endpoint on the stretched horizon : 



1 



T**(z)d z ll>r{u,z) 



z=z h 



(3.14) 



We now analyze the consequences of these equations and, in particular, emphasize the 
differences w.r.t. the corresponding analysis in Ref. S. 



Eq. (3.14) is a Langevin equation of a special type: the noise term is precisely compen- 
sating the pulling force T zz (zf l )d z ip r due to the string tension. By taking the expectation 
value of this equation and recalling that T zz (zh) ~ e vanishes in the limit e — > 0, we 
conclude that d z (i(^(uj, z)) is regular near the world-sheet horizon; this implies that the 
average value of the classical solution is proportional to ipret : 



(ip r (u,z)) = (i$(oj))il> re t(uj,z) 



(3.15) 

the boundary value of 



The normalization is fixed by the expectation value of ^{oo) 
tp r (u,z) at z = z m <C 1. 

We will now construct the solution ip r {u), z) to the EOM ( |2.7| ) by specifying its bound- 
ary values, ip^(uj) and ip^(u>), at the endpoints z = z m and z = z^, respectively. After 
simple algebra, the respective solution can be written as 



i/j ret (uj,z) - j> adv (uj,z) 
i/; ret (uj,z h ) - ip adv (u,z h ) 



(3.16) 



The reason why this particular writing is natural is as follows: when taking the expectation 



value according to Eq. ( 3.15 ), we find {^{ui)) = (tp^(uj)) ip ret (uj, Zh), which shows that the 
coefficient ^i 00 ) ~ V'?( w )V ; ret(w ) Zh) in front of the second term in Eq. fl3,16|) is a random 
variable with zero expectation value. Clearly, this term plays the role of a noise. The 
statistics of this noise is determined by the horizon Langevin equation ( |3.14j ), and in turn 
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it implies a boundary Langevin equation for i/j®(u>), via the condition ( |3,13 ), Let's see how 
all that works in detail. We will first rewrite Eq. ( |3.16[ ) as 

1p r (uJ,z) = i%(uj)ll) re t(uj,z) +lC(U0) — -, r , (3.17) 

thus fixing the normalization of the noise term £°(w). After inserting Eq. ( |3.17| ) in the 
Neumann boundary condition ( |3.13| ) (say, in the transverse sector), one finds 4 

G° R (u)^ r (u;) =eV), (3-18) 



which is the standard form of a Langevin equation (compare to Eq. (|2,36| )). 

It remains to check that the statistics of £°, as inferred from Eq. ( |3.14| ), is indeed the 
same as previously derived in Sect. 2.3. To that aim, we insert the form (3.17) of the 
solution into Eq. ( p. 14 ); as already explained, the regular piece of the solution oc ifi re t(u}, z) 
does not contribute in the limit e — > 0, so we are left with 

- 5 d^ r " (z) **- (a, ' z) U = <*(«)■ < 3 - 19 > 

\JsingT zz (z)d z ip adv (uj, z)\ z=Zh = -2iujjr]ip adv (uj, z h ), d. Eq. (|3j|) , together with Im Gr(u) = 
—ujj7], this finally becomes 

^adv(u},z h )^(u) = £ h (u). (3.20) 

This relation is in fact natural, as we argue now: from the transformation connecting £ 
to if) a , or directly by comparing Eq. ( p. 17 ) with the standard expression (2.22) for rf) ri one 
can see that the strength of the noise term scales like £(cj) ~ ip a (w) G sym (u>). On the other 
hand, Eq. ( f2.22| ) implies ipa(uj) = ip adv (uj, Zh)ip a {uj)- Hence one can write 

£ h H ~ VaMG s h ym M - lpadv(u,Z h )^ a (L0)G syra {L0) ~ ^ adv (u, Z h ) £° (w) . (3.21) 

Remarkably, the relative factor between £° and £ h in Eq. ( |3.2C| ) does not spoil the nor- 
malization of the noise-noise correlator, because \ipadv{^i z h)\ = 1 ; as we now demon- 
strate. To that aim we rely on the observation at the end of Sect. 2.2 that the r.h.s. of 
Eq. ( |2.32| ) is independent of z. Clearly, this remains true after replacing ip ret — ► ipadv i n 
Eq. ( |2.32; ). (Indeed, the current ( p. 33 ) is conserved for an arbitrary solution ^f so i(t,z).) 
Writing ip adv (uj, z) = C(uj)^ dv (uj, z), so that ip adv (uj, z^) = C(u), and evaluating the r.h.s. 
of Eq. ( |2.32|) separately at z = z m and z = Zh, one deduces that |C(o;)| = 1, as anticipated. 
Thus, the 2-point function (£°(u;)£ (u/)) is indeed the same as in Sect. 2.3. Note also that 
our previous argument is independent of the precise value of e (the distance between the 
world-sheet and the stretched horizons), so long as e is small enough for the near-horizon 
expansions to make sense. This strongly suggests that the strength of the noise remains 
constant along the string, from the stretched horizon up to the boundary. 



4 



The following identities, which can be checked from Eq. (2.25), are useful in this respect: 

Gu,ij(w) = -^T±(Zm) [d z lp r£t {u,z) + dzlpadv{-U,z)\ z = Zm . 

ImGfl(oj) = - T±(z m ) d z [V>ret(a>, z) - ip adv (uj, z)] z=z ^. 
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4. Discussion and physical picture 



In this section, we will first discuss some consequences of the previously derived Langevin 
equations, which support the idea that the noise terms in these equations are of non- 
thermal nature, and then propose a physical picture in which these fluctuations are inter- 
preted as quantum mechanical fluctuations associated with medium-induced radiation. 

4.1 Momentum distributions from the Langevin equations 

An important property of the Langevin equations (^il )-(|2.45|) that we would like to em- 
phasize is that, except in the non-relativistic limit 7 — 1, these equations do not describe 
the thermalization of the heavy quark. There are several arguments to support this conclu- 
sion. For instance, in thermal equilibrium the momentum distributions should be isotropic, 
but this is clearly not the case for the large-time distributions generated by Eqs. ( |2.44 )- 



(2.45), because of the mismatch between K£ and k± when 7 > 1. Besides, in order to 



generate the canonical distribution for a relativistic particle, P oc exp{ — yp 2 + Mq/T}, 
the noise correlations must not only be isotropic, but also obey the relativistic version of 
the Einstein relation, which reads k = 2ETrj£, [ pq| . Using Eqs. (2.37) and ( p. 46 ), it is easily 



seen that this condition is not satisfied for either transverse, or longitudinal, fluctuations 
(except if 7 = 1, once again). 

The Einstein relation is a particular form of the fluctuation-dissipation theorem, so its 
failure might look surprising given that the Green's functions at the basis of our Langevin 
equations obey the KMS condition ( p. 26 ). Recall, however, that this peculiar KMS condi- 



tion involves an effective temperature T e fj = T/7 1 / 2 ; and indeed, in the transverse sector 
at least, the Einstein relation appears to be formally satisfied, but with T — > T e ^. But this 
does not hold in the longitudinal sector, where K£ involves an additional factor 7 2 . Hence, 
the present equations cannot lead to thermal distributions. 

It is then interesting to compute the actual momentum distributions generated by these 
Langevin equations at large times. Consider first the transverse sector, and introduce the 
probability distribution P(p±,t) for the transverse momentum p± = (pi,P2) at time t: 

p(p±,t) = f m\5(p L - Pl _u{t)) e~^ Atmm . (4.1) 



Here, p±[Ci]{t) is the solution to Eq. ( [2.44 ) corresponding to a given realization of the 
noise 5 , and reads (with i = 1, 2; we assume Pi(0) = so that (p±(t)) = at any time) 

Pl (t) = fdt'e-^^Ut')- (4.2) 
This implies (p\) = (p 2 ,) = with 







<£(*)> = - ^, (4.3) 

2rj D 2r] D 



In general, this solution will depend on our prescription for discretizing the time axis; this is so since 
the noise-noise correlator depends itself on the momentum ('multiplicative noise'). But to the accuracy 



of interest, we can treat 7 in Eq. (2.37) as the fixed quantity 70, and then one can safely use continuous 
notations. 
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where the last, approximate equality holds for large times r]Y)t S> 1. Returning to Eq. ([4.11), 
this gives 

^-ssbjH-^i}- (44) 

A similar expression holds in the longitudinal sector, but only after subtracting away the 
global motion with velocity Vq, which one can do by writing 5pg = P£—po with po = Mq'JqVq. 

For large times t 3> 1/?7d> the transverse and longitudinal momentum distributions for 
the heavy quark approach the following, stationary, forms 

1 / pi + p\ 



P(p±,t) ~ - — 1/9 __. exp 



2vr7 1 /2rM Q I 2j 1 / 2 TMq 



1 



P(6l)l ./) -~ - exp{- | , (4.5) 



which formally look like thermal, Maxwell-Boltzmann, distributions for non-relativistic 
particles, but with different temperatures in the transverse and longitudinal sector - 
Tj_ = 7 1 / 2 T and Tg = 7 5 / 2 T — , none of them equal to the plasma temperature T. 

4.2 Physical picture: Medium induced radiation 

In this section, we propose a physical picture for the dynamics of the heavy quark, as 
encoded in the Langevin equations ( 2.44 )— fl2.46|) . The general picture is that of medium- 



induced parton branching, as previously developed in Refs. (12], [l3|, |l4|, [|, that we shall 
here adapt to the problem at hand. As we will see, this qualitative and admittedly crude 
picture provides the right parametric estimates for both the drag force and the (transverse 
and longitudinal) momentum broadening. Besides, it supports the non-thermal nature of 
the noise terms in the Langevin equations. 

Due to its interactions with the strongly-coupled plasma, a heavy quark can radiate 
massless N = 4 SYM quanta (gluons, adjoint scalars and fermions) which then escape in 
the medium, thus entailing energy loss towards the plasma and momentum broadening 
(due to the recoil of the heavy quark associated with successive parton emissions). This 
dynamics is illustrated in Fig. |2[ The main ingredients underlying our physical picture are 
as follows: 

(i) The emission of a virtual parton with energy u> and (space-like) virtuality Q 2 = 
k 2 — ui 2 > requires a formation time i co h ~ u>/Q 2 . (k is the parton 3-momentum, and 
we assume high-energy kinematics: \k\ ~ u 3> Q-) This follows from the uncertainty 
principle: in a comoving frame where the parton has zero momentum, its formation time 
is of order l/Q; this becomes uj/Q 2 after boosting by the parton Lorentz factor 7 P = u/Q. 
Note also that, when the parent heavy quark is highly energetic (7 3> 1), the momentum k 
of the emitted parton is predominantly longitudinal: |fc| ~ kg ~ uj, whereas k± ~ Q <S kg,. 

(ii) During the formation time t co ^ ~ ui/Q 2 , the heavy quark does not radiate just 
a single parton, but rather a large number of quanta, of 0(y\), whose emissions are 
uncor related with each other. This is merely an assumption, which as we shall see provides 
the right A-dependence for the final results. 
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(iii) Only those quanta can be lost towards the plasma, whose virtualities are small 
enough — smaller than the saturation momentum Q s ~ t co hT 2 corresponding to the parton 
formation time. This follows from the analysis in Ref. [^] which shows that an energetic 
parton propagating through the strongly coupled plasma feels the latter as a constant force 
~ T 2 which acts towards reducing its transverse momentum (or virtuality). Then, a space- 
like parton, which would be stable in the vacuum, can decay inside the plasma provided 
the lifetime t co h of its virtual fluctuations is large enough for the mechanical work ~ t co hT 2 
done by the plasma to compensate the energy deficit ~ Q of the parton. This condition 
amounts to Q < Q s , with the upper limit given by 

Qs ~ t coh T 2 ~ (u/T 2 ) 1 / 3 ~ ^%T, (4.6) 

where we have also used t co h ~ uj/Q 2 and 7 P = uj/Q. 

(iv) The rapidities of the radiated quanta are bounded by the rapidity of the heavy 
quark: j p < 7. This is again motivated by the uncertainty principle and at least at weak 
coupling it is confirmed by the explicit construction of the heavy quark wavefunction [^6| . 

We shall now use this picture to compute the rate for energy loss and momentum 
broadening of the heavy quark. The latter radiates energy AE ~ y/Xoj over a time interval 
At ~ uj/Q 2 , where uj and Q are constrained by Q < Q s (u>,T). The dominant contribution 
to the rate \AE/At\ comes from those quanta carrying the maximal possible energy uj ~ "fQ 
and also the maximal corresponding virtuality Q ~ Qs{l, T) ~ \flT (to minimize the 
emission time). Therefore, 



dE x/X 



uJ 



VXQ 2 ~y/X 7 T 2 , (4.7) 



dt (u/Q 2 ) 

in qualitative agreement with the estimate for the drag force -Fdrag = r ll v ~ 7V^XT 2 in 



Eq. ( 2.45 ). (Recall that we consider the relativistic case v ~ 1.) 

Consider similarly momentum broadening: being uncorrelated with each other, the y/\ 
quanta emitted during a time interval £ co h have transverse momenta which are randomly 
oriented, so their emission cannot change the average transverse momentum of the heavy 
quark. However, the changes in the squared momentum add incoherently with each other, 
thus yielding (once again, the dominant contribution comes from quanta with Q ~ Q s (l, T) 
and uj ~ 7Q) 

d(pi> VXQ 2 r-Qi 1— — o 

~ ; u 7 m2\ ~ vX— VX^fyT , 4.8) 

dt (u/Qi) iQs 

which is parametrically the same as the estimate for k± in the first equation fl2.37|) . The 
random emissions also introduce fluctuations in the energy (or longitudinal momentum) of 
the heavy quark, in addition to the average energy loss. The dispersion associated with 
such fluctuations is estimated similarly to Eq. (|4.8|) (below, 5pg = pi — (pi)) 



dt (w/Qi) 



-21- 



Figure 2: Energy loss and momentum broadening via medium-induced parton emission at strong 
coupling. It is understood that the radiated partons feel a plasma force which allows them to be 
liberated from the parent heavy quark (see text for details). 



in qualitative agreement with the previous result, Eq. ( 2.3?1 ), for k,£. Note that, with this 
interpretation, the relative factor 7 2 in between ki and k± is simply the consequence of 
the relation oj ~ jQ between the energy and the virtuality (or transverse momentum) of 
an emitted parton. 

This physical picture also clarifies the role of the world-sheet horizon in the dual gravity 
calculation: via the UV/IR correspondence, the radial position z s = 1/^/7 of this horizon 
(in units of zh = V 71 "^) * s mapped onto the saturation momentum Q s ~ ^yT in the 
boundary, gauge, theory. Hence the emergence of the noise terms from the near-horizon 
dynamics of the string reflects quantum-mechanical fluctuations in the emission of quanta 
with virtualities Q ~ Q s , which as we have just seen control momentum broadening. 

It is furthermore interesting to compare the above physical picture to the corresponding 
one at weak coupling [|l5| , |l6| , [T^ , |i"8| . Note first that the mechanism for momentum 
broadening is different in the two cases: at weak coupling, this is dominated by thermal 
rescattering, i.e., by successive collisions with the plasma constituents which are thermally 
distributed (see Fig. |3[). In that case, the rate d(p^_)/dt = q defines a genuine transport 
coefficient — the "jet-quenching parameter" — , i.e. a local quantity which depends only 
upon the local density of thermal constituents (quarks and gluons) together with the gluon 
distribution produced via their high-energy evolution. By contrast, at strong coupling, 
the dominant mechanism at work is medium-induced radiation, which is intrinsically non- 
local (it requires the formation time t C oh) an d hence cannot be expressed in terms of a local 
transport coefficient. Medium-induced radiation is of course possible at weak coupling too 
(see Fig. |j), but the respective contribution is suppressed by a factor g 2 N c as compared to 
the thermal rescattering. We see that, formally, it is the replacement g 2 N c — ► y/\ (i.e., the 
coherent emission of a large number of quanta) which makes the medium-induced radiation 
become the dominant mechanism for momentum broadening at strong coupling. 

On the other hand, energy loss is predominantly due to medium-induced radiation 
at both weak and strong coupling, but important differences occur between the detailed 
mechanisms in the two cases (compare Figs. |2| and ||): At weak coupling, the radiated 
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Figure 3: Momentum broadening via thermal rescattering at weak coupling. 




Figure 4: Energy loss via medium-induced gluon emission at weak coupling 



gluon, which typically comes from a highly virtual gluon in the quark wavefunction, is 
freed (radiated) via thermal rescattering. At strong coupling, radiation is caused by the 
plasma force ~ T 2 . After being emitted, the parton undergoes successive medium-induced 
branchings, thus producing a system of partons with lower and lower energies and transverse 
momenta, down to values of 0(T), when the partons cannot be distinguished anymore from 
the thermal bath. 

It is finally interesting to notice that, in spite of such physical dissymmetry, the formula 



for energy loss at weak coupling can be written in a form which ressembles Eq. (4.7), namely 

dE 

— ~ g 2 N c Q 2 s (weak coupling) , (4.10) 

where however Q s is now the saturation momentum to lowest order in perturbative QCD 
and is related to the respective jet-quenching parameter via Q 2 ~ t/i co h- Energy loss 
involves a coherent phenomenon at both weak and strong coupling. 
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